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Abstract 

We review the Hall-Vinen-Bekarevich-Khalatnikov (HVBK) equations for superfluid 
Helium turbulence and discuss their implications for recent measurements of superfluid 
turbulence decay. 

A new Hamiltonian formulation of these equations renormalizes the vortex line 
velocity to incorporate finite temperature effects. These effects also renormalize the 
coupling constant in the mutual friction force between the superfluid and normal fluid 
components by a factor of ps/ p (the superfluid mass fraction) but they leave the vortex 
line tension unaffected. Thus, the original HVBK form is recovered at zero temperature 
and its mutual friction coefficients are renormalized at nonzero temperature. The 
HVBK equations keep their form and no new parameters are added. However, a 
temperature dependent trade-off does arise between the mutual friction coupling and 
the vortex line tension. 

The renormalized HVBK equations obtained via this new Hamiltonian approach 
imply a dynamical equation for the space-integrated vortex tangle length, which is the 
quantity measured by second sound attenuation experiments in superfluid turbulence. 
A Taylor-Proudman theorem also emerges for the superfluid vortices that shows the 
steady vortex line velocity becomes columnar under rapid rotation. 

PACS numbers: 47. 10. -Kg, 45.10.Db 
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1 HVBK equations 

Recent experiments establish the Hall-Vinen-Bekarevich-Khalatnikov (HVBK) equations as 
a leading model for describing superfluid Helium turbulence. See Nemirovskii and Fiszdon 
[1995] and Donnelly [1999] for authoritative reviews. See Henderson and Barenghi [2000] for 
a recent fluid mechanics study of steady cylindrical Couette flow using computer simulations 
of the incompressible HVBK equations. 

In the Galilean frame of the normal fluid with velocity v„, the HVBK equations may be 
expressed as follows, upon ignoring thermal diffusivity and viscosity, 

dtp = - divJ , dtJ, = -dk T> , dtS=- div(Sv„) + R/T , 

PsdtVs + psi-Vs ■ V)Vs = -p^V(yU - -|V^ - V„P) + p^f , (1) 

and summing over pair of upper and lower repeated indices. One may consult, e.g., Bekare- 
vich and Khalatnikov (BK) [1961] and Donnelly [1999] to compare these equations with the 
form they take in the reference frame of the superfluid.^ 

Notation. Here p and ps denote the total and superfluid mass densities, respectively. The 
entropy density of the normal fluid is S, its temperature is denoted T, and Pn = p — Ps 
denotes its mass density. The superfluid velocity is denoted and 

J = Ps^s + Pn^n 

is the total momentum density. In the entropy equation R is the rate that heat is produced 
by the phenomenological friction and reactive forces in f , which must be specified to close 
the theory. We also denote 



Stress tensor 
Euler's pressure law 
Superfluid First Law 



T^^ = Ps VsiV^s + Pn Vni + (P + A-Cj) 5f - XiU^ , 

P = -Eq + TS + pp, 

dsQ = pdp + TdS + ( J — pv„) ■ d{-Vs — v„) + A ■ (ia; . (2) 



In the superfluid First Law, uj = curl Vg is the superfluid vorticity with magnitude \uj\ = Cj-uj, 
where Cj = uj/\uj\ is its unit vector. BK [1961] takes the energy density Eq to depend on the 
magnitude of the superfluid vorticity, Ic^l, as 

£0 = — In - . (3) 

47r a 

This is the energy per unit length of a superfluid vortex line, ps(K^/47r)ln(/2/a), with quantum 
of circulation k = h/m ~ 10~^(cm^/sec) and ratio R/a of mean distance between vortices R 
to effective vortex radius a, times the vortex length per unit volume, \uj\/k. Hence, we flnd 

der, d\uj\ , , , „ . , , , , , K ^ R 

^ = 7^^ = ^^^ A/p, ^— ln- = Ao, 4 

d\uj\ ouj A-K a 



^In making this comparison it is useful to recall the Galilean transformation of the chemical potential, 
/i' = /i — i|vs — v„p , where fi' is evaluated in the superfluid frame and /i in the normal-fluid frame. See 
Puttcrman [1974] for a clear discussion of the role of Galilean transformations in superfluid hydrodynamics. 
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where one ignores the derivative of R ~ ^Jk/\uj\ inside the logarithm. The appearance of 
A in the stress tensor Tj' shifts the pressure P, and divT introduces an additional force 
— uj ■ VA = — PsT into the motion equation. The quantity T is called the "vortex line 
tension." 

BK [1961] assigned the following form to the phenomenological coupling force f appearing 
in the superfluid velocity equation in (|1|), 

f = (vl-v^)xu;, Vi = v„ - (as° + /5c2; X s°) , (5) 
where s° = v° — v„ , v° = Vs + p^^'^curl A , (6) 

with "vortex velocity" and "slip velocity" s° introduced as auxiliary quantities. The 
HVBK equations written as (|l|) in the normal-fluid frame conserve the energy, 



E 



"""plvnl^ + (J - pv„)-v„ + eo]c?^a;. (7) 



-2 

The BK [1961] form of the phenomenological force f also implies the dissipative heating rate, 

R= {3 - pwn + curl A) ■ X (vl - v„) , (8) 

which is Galilean invariant and positive. Substituting this form of f into the superfluid 
motion equation in (|I|) and taking its curl provides the following equation for the superfluid 
vorticity, u = curlvg, 

dtio = curl (vl X uo) . (9) 
This vorticity equation implies the HVBK superfluid Kelvin circulation theorem, 

^ I ( u-dS = ^i v,-dx = 0. (10) 
dt J Js dt Jds{vL) 

The Kelvin formula (|T0|) expresses conservation of the flux of superfluid vorticity through 
any surface S whose boundary dS moves with the velocity v^, so may be regarded as 
the local velocity of a vortex line. Equivalently, the Kelvin formula expresses conservation 
of superfluid velocity circulation around any loop that moves with the vortex line velocity 

In a key phenomenological step that closed the theory, BK [1961] assigned the undeter- 
mined functions a and (3 in the force f and its auxiliary vortex line velocity as 

l + ap. = ^, PPs = ^. (11) 

The dimensionless coefficients B and B' were introduced earlier in Hall and Vinen (HV) 
[1956] to parameterize the Gorter-Mellink [1949] "mutual friction" force (B) and its reactive 
component {B'). Hence the name, HVBK equations for this closure. 

The assignments in BK [1961] of the undetermined functions a, [3 in ([Tl|) , as well as 
the vortex slip velocity s° in (^) are designed to reproduce the phenomena observed in HV 



Renormalized HVBK dynamics for He-II 



D. D. Holm 



4 



[1956] and yet still conserve mass, momentum and energy. This phenomenological approach 
used in BK [1961] is indeterminate, however, in the sense that some freedom still remains 
in making these assignments. The HVBK equations (|l|) that result from this approach 
do possess the desired conservation laws for mass, momentum and energy. And they also 
possess a Kelvin theorem for the circulation of superfluid velocity. However, because of the 
indeterminacy inherent in the phenomenological approach, the HVBK equations (|l]) are not 
unique in possessing these properties. An alternative assignment of the vortex slip velocity 
is 

s = V£ — v„ with = V + p^^curl A , with v = p~^(psVs + p„v„) . (12) 

The mean velocity v also figures prominantly in Hills and Roberts [1977] discussion of the 
HVBK equations. As we shall see, the alternative expression ( [T^ ) for the auxiliary vortex 
slip velocity in terms of the mean velocity v arises naturally in the Hamiltonian derivation 
of a slightly modified set of HVBK equations. These equations possess the same formal 
conservation and circulation properties as HVBK, modulo redefining the vortex slip velocity 
as s rather than s°. The vortex slip velocity s in ( [T^ ) is defined relative to the Galilean frame 
of the normal fiuid, which is present only at finite temperature. The HVBK s" in (^ is the 
limit of the vortex slip velocity s for zero temperature, at which no normal fiuid remains. 
The vortex slip velocity s in (|T2D is a slight modification of s° in (^ necessary to incorporate 
finite temperature effects, without changing the form of the HVBK theory, obtained from 
a Hamiltonian derivation of these equations in the normal fluid frame. In the Hamiltonian 
framework, the energy-momentum conservation laws and Kelvin circulation theorem are all 
natual consequences. Moreover, the velocities and v„ are identifled as being dual to the 
momenta given by pv^ and Pn(v„ — v^), respectively. 

Outline. We shall use a Hamiltonian approach with Lie-Poisson brackets to derive the 
expression (p!^ ) for the vortex line velocity at flnite temperature from flrst principles by 
using the energy E in (^ as the Hamiltonian. The momenta conjugate to the velocities 
and v„ shall be our basic dynamical variables. The flnite temperature vortex line velocity 
Ye and slip velocity s determined this way turn out to be 

= p"^ (J + curl A) , and s = - v„ =i (— ) (v^ - v„) . (13) 

V p ^ 

At zero temperature, p —>■ ps and these reduce to the BK [1961] phenomenological expressions 
with v° given by (^. Thus, the flnite temperature corrections found by using the Hamiltonian 
approach renormalize the HVBK slip velocity in the mutual friction force f by the factor 
ps/p (the superfluid mass fraction). Aside from this renormalization, the vortex line tension 
is left unaffected by this renormalization, the superfluid vortex equation keeps its form and 
no new parameters are added. 

Technical details of deriving this renormalized theory from its Hamiltonian and Lie- 
Poisson brackets are given in the Appendix. 
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Main results. We shall use the superfluid vortex dynamics for the renormahzed HVBK 
equations obtained via this Hamiltonian approach to write a dynamical equation for the 
space-integrated total vortex tangle length, which is the quantity measured in the Oregon 
experiments on superfluid turbulence reported in Skrbek, Niemela and Donnelly [1999]. 

We shall also study the restriction of the renormahzed HVBK equations for the incom- 
pressible case, in which p„ and ps are constants and one takes V • v„ = and V • Vs = 0. 
Finally, we shall demonstrate the invariance of the forms of these equations upon transform- 
ing into a rotating frame. The Coriolis force in such a rotating frame couples to the vortex 
line velocity v^, which of course differs from both the superfluid velocity and the normal 
velocity. We shall derive a Taylor-Proudman theorem for steady superfluid vortices un- 
der rapid rotation. According to this superfluid Taylor-Proudman theorem, the vortex line 
velocity becomes columnar under sufficiently rapid rotation. That is, the lateral vortex line 
velocity is nondivergent and independent of the axial coordinate, and the axial velocity de- 
couples from the lateral motion. Therefore, under sufficiently rapid rotation, the superffuid 
vortex filaments will straighten and become parallel to the axis of rotation as they approach 
a steady state. 

Numerical implications. This renormalization of the vortex line element slip velocity 
in the HVBK equations from s° — ^ s ~ s^ps/p is sensitive to temperature, but it does not 
affect the vortex line tension. Therefore, a temperature sensitive trade-off arises between 
mutual friction and vortex line tension that may be worth testing in numerical simulations 
such as those reported in Henderson and Barenghi [2000] . The HVBK equations are thought 
to break down in the presence of strong counterflow. However, as general conservation 
laws there is no mechanism in the equations that would signal this breakdown. A rotating 
Rayleigh-Besnard experiment might be useful in testing the range of validity of the HVBK 
equations (Barenghi, private communication). Such an experiment might also indicate how 
these equations should be modifled in the presence of strong counterflow. 

Experimental implications. Temperature sensitivity of the coupling between the super- 
fluid vortices and the normal fluid component is an area of intense current investigation in 
superfluid Helium turbulence, see Donnelly [1999]. One would like to know whether the 
Ps/p renormalization of the mutual friction forces relative to the vortex line tension would 
matter signiflcantly in comparisons of the predictions of the HVBK equations with modern 
experiments in Helium turbulence at low, but flnite temperatures. 

Superfluid vortex dynamics. To begin addressing this issue, we may use the superfluid 
vorticity equation for the renormalized HVBK equations obtained in the Appendix via the 
Hamiltonian approach to write an explicit equation for the dynamics of Vinen's vortex 
length density L = \lv\/k. In the superfluid turbulence decay experiments reported by 
Skrbek, Niemela and Donnelly [1999] the spatial integral of this quantity is measured as a 
function of time to decrease over six decades as The integrated vortex length measured 
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in these experiments is predicted by the renormahzed HVBK equations to be governed by 
the superfluid vorticity dynamics alone. 

Upon including mutual friction, the superfluid vortex dynamics for the renormalized 
HVBK equations is expressed as, cf. equation @, 

dfUJ = curl(vi X uj) , (14) 

in which the renormalized total vortex line velocity given by, cf. equation (^, 

B'pn Bpn 

vl = ^£ —uxs, where s = V£ - v„ , (15) 

and its Hamiltonian limit is found to be 

= V + p^-'^curl A , with v = p^^J and X = \uj. (16) 

Thus, relative to the Hamiltonian approach, the terms in B and B' are additional velocities 
introduced by phenomenology, while V£ is the vortex line velocity in the absence of mutual 
friction. 

The HVBK superfluid vorticity equation implies the following dynamics for the integrated 
vortex length measured in the turbulence decay experiments, 

UJ ■ dtuj / K d^x 




vorticity dynamics 

= J Lv ■ (6j X curlii) d^x — J L {uj x cuyW uj) ■ {uj x curio)) d^x 

transport ■ curvature damping by curvature 

+ <f L{nxij)- [{Ljxv,) + ^{ve~v,))dS . (17) 
^ ^ 
creation and destruction at the boundary 

Here uj is the unit vector tangent to a superfluid vortex filament, so k, = {uj x curlu)) is its 
local curvature. The transport and damping of the vortex tangle length is proportional to 
this local curvature. The effective velocity v in the transport term is given by 

_ B'pn . Bpn . ,_ . , . 

v = v-^— (v-v„)-^o;x (v-v„). (18) 

According to the last term in (|l^, vortex length is created or destroyed at the boundary, 
unless the vortex filaments approach it in the normal direction, so that fi x cl; = 0. 

Formula (0) for the evolution of the total superfluid vortex length presents a trade-off 
between the mass- weighted velocity v and the local induction velocity (or filament curvature) 
UJ X curio), in the interior of the domain. This trade-off in the interior competes with the 
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process of creation and destruction at the boundary. For example, in counterflow turbulence, 
the superfluid moves toward the heater at the boundary, so the term in v would tend to be 
nonzero. In contrast, for grid turbulence, v is small, so this term would tend to contribute 
less. This formula governs the dynamics of the experimentally measured quantity / Ld^x. 
However, it does not yet show how to obtain the decrease seen in this quantity by 

Skrbek, Niemela and Donnelly [1999] in their experiments on decay of turbulence. 

Suppose the main source of decay were the term labeled "damping by curvature" in 
formula (0) and the flow were isothermal and incompressible. This would imply 

1 V>--o(r)A„<^/;^!l^-^. (19) 



{L)dV ' ' " {L) 2(t + to)' 

where Aq = X/ps = / 4:Ti)ln{h / a) is the quantum vortex constant, to is a time shift in the 
experimental analysis, Cq{T) = BpnPs/{2p'^) and angle brackets (■) denote spatial integral 
over the measurement domain. In particular, 

(L) = J Ld^x, {L/R^) = J L\ujx cnr\uj\'^d^x. (20) 



The measured decrease in (L) implies via formula ( |19|) that the length-weighted mean 
curvature of the vortex tangle {L/R^)/{L) decays due to mutual friction as . Thus, on 
the average as the vortex length decays, the vortices tend to straighten, under the 
effects of mutual friction damping. 



Preservation of helicity versus preservation of vortex length. The helicity, or 
linkage number for the superfluid vorticity is defined as 

A = j{ws-uj) d^x. (21) 

The helicity satisfies an evolution equation obtained from the superfluid vortex dynamics, 

— = - j> {h ■ uj)(^p - -vl - Vs ■ {vl - v„)) dS - j>{n- Vl) (v^ ■ uj) dS . (22) 

Therefore, even with mutual friction, helicity is created and destroyed only on the boundary. 
Moreover, helicity will be preserved, provided both u and v^ are tangential at the bound- 
ary. The former condition, however, is the opposite of that required for the creation and 
destruction of vortex length at the boundary to cease. Therefore, no equilibrium should be 
expected that preserves both the helicity and the vortex length in a superfluid. 



Superfluid vortex equilibria are not ABC flows. The steady equilibrium solutions of 
the superfluid vorticity dynamics satisfy 

curl (vl X cj) = . (23) 

For example, a steady equilibrium exists when uj and v^ are parallel. Note that these 
"super-Beltrami flows" are not eigenf unctions of the curl. Therefore, they are not Arnold- 
Beltrami-Childress (ABC) flows, as occur for the Euler equations. 
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2 Incompressible renormalized HVBK flows 

To express the renormalized HVBK equations in the incompressible limit, we begin by rec- 
ollecting the compressible equations and abbreviating jv^ — v„p = v^^ , 

dtS = -div(5v„,) + i?/T, 
dtp = - div(p^v^ + p„v„) , 

ps(5tv, + (v, ■ V)vJ = -p^V(/i- +Ps(vl- V,) X cu, 

dt{PsVsi + PnVni) = - <9fc {PnVniVn + PsVsiv'l) - OiP - dkT^ , 

As we have seen, finite temperature effects renormalize the total vortex line velocity as 

Vl = — — ( — u) X s H , where s = — v„ , (24) 



p\2 2 

and the Hamiltonian part of the line velocity (with corrections for finite temperature) is 
defined as 

V£ = V + p^-'^curl A , with v = p^^(psVs + p„v„) and A = Aa;. (25) 

To the extent that p, p^, p„ and S all may be taken as constants for a given temperature 
and the heating rate R is negligible, then the velocities v„ and are incompressible, i.e., 

V ■ v„ = and V ■ v, = . (26) 

In this situation, the pressure P may be obtained from the Poisson equation, 

- V'(P + \-uj) = div(p,(v, ■ V) V, + p„,(v„ • V) v„, - ■ VA) , (27) 

found from the divergence of the total momentum equation. Combining the superfiuid 
motion equation with total momentum conservation results in an equation for the normal 
fluid velocity in the incompressible case 

Pndt^n + Pn(v„ ■ V)v„ = - V( P' - PsP + ^t^L ) " Ps (vl - V,) X ^ + ^ ■ VA , (28) 



where is given in equation (^4]). We set P' = P + A ■ and take it as the total pressure. 
(One also could have absorbed A ■ uj into P earlier, by including it in Euler's pressure law.) 
Since A = |A|a; and a; is a unit vector, we find for constant ps the standard relation for the 
vortex line tension denoted as T. Namely, 

uj ■ VA = — Ao Ps X curl Cj = psT , (29) 

where Aq = X/ps = (K/47r)ln(6/a) is a constant. 
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Remark. We note that the quantity T known as the vortex hne tension first appears in 
the normal Buid equation, as a reaction to the presence of the superfluid. The standard 
convention for introducing the mutual friction force has the effect of shifting T into the 
superfluid equation. By action and reaction, though, T could appear in either equation. 

These equations of motion must be completed by providing an equation of state relation 
for the quantity p — \vl^. BK [1961] assumes a law of partial pressures, 

P^=P^P' = P'-P, and P, = ^P' = p,^-lp,T;2^. (30) 
P P 2 

In this case, the renormalized HVBK motion equations for incompressible flow reduce to 

9tv, + (v, ■ V)v, = -lvP'-^F„, + T, (31) 

P P 

5tv„ + (v„ ■ V)v„ = -ivP' + ^F„,. (32) 

In these superfluid motion equations, the renormalized mutual friction force F^,; is 
defined as the sum (with uj = curlvs) 

F„, = (s°x.;) + (^) FO,, where F", = (| .2; x s° + ^ s°) x ^ . (33) 

Here F°^ is the HVBK mutual friction force without any finite temperature corrections. To 
acquire these formulas, we used the relations for the incompressible case, 

s = V£ - v„ = — (vs + Ao curia; - v„) = — (v° - v„) = — s° , (34) 
P P P 

with s'' = Vs + Ao cutIlj — v„ , and we eliminated by using the relation 

- Ps(vl - - AoCurlLj) = Pn s + ^-^(y ^ ^ ^ "Y ^) ■ ^^^^ 

In equation (^) for F„s, the quantity {s^ x u) is the Hamiltonian reactive force (which 
could be naturally absorbed into Vinen's B' parameter) and F°^ is the phenomenological 
mutual friction force defined according to the standard convention as in BK [1961] and 
Donnelly [1999]. The finite-temperature corrections contribute an overall factor of Ps/p 
to the standard zero-temperature expression F^^ for the phenomenological mutual friction 
force. No new parameters are added, but a temperature dependent trade-off is identified 
between the renormalized mutual friction coupling and the vortex line tension, since the 
vortex line tension remains unaffected by the finite-temperature corrections. 

In the isothermal case, the motion equations are closed by the Poisson equation for P, 
since the other coefficients (P, B', Pn/ p, etc.) are specified functions of temperature and 
they may be taken as constants, for an isothermal incompressible superfluid flow. 
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Note that equations (|3Tl - p2D may be rewritten with Un = curl v„ as 

dt^^s + Vfis = v„xcu + ^s°xa;-^F°, with p, = P'/p + It;^ , (36) 

atV„ + V/i„ = v„xcu„ + ^s°xcu + %Fl with fi^ = P'/p+^vl. (37) 

P P^ 2 

These equations imply an equation for the velocity difference, 

5t(v, - v„) + lv{vl - vl) = V, X (cu - cu„) - ^F°, with cu, = curl . (38) 

and there is no tendency for mutual friction to cause any alignment in the vorticities of 
the superfluid and its normal component. Instead, the curlcj part of F°^ 7^ would break 
any such alignment, if it were to form spontaneously. Indeed, alignments sufficient for 
steady solutions are 

s° X (V/is X V/i„) = , s° X = and s° X v„, = , with s° = v^ + Aocurlcj — v„ , (39) 

provided ps and p„ are functionally unrelated. Thus, the steady equilibrium alignments 
imposed by mutual friction involve v^, and curlci), as well as the independent gradients 
of ps and pn. For example, one class of equilibria has s°, v„, lo all aligned tangent to 
intersections of level surfaces of pn and ps- 

3 Rotating frame renormalized HVBK equations 

We transform to a rotating frame with relative velocities denoted with an asterisk as v* = 
v^ — R(x), etc., and curlR = %\. After a calculation involving Legendre transformations, we 
obtain the Hamiltonian for the relative motion, cf. the Hamiltonian in (plf ) of the Appendix, 

MM*,p,^,u,A*,n) = ||lp|v: + R(x)|2 + (M*-pA*-pv:)-v: (40) 
+ £o(p, v: - v:, 0^* + 2fi) - R ■ [p(v: + R) + (p - n)(A* - R)] jd^^; . 

Here M* — pA* = J* = J — pR and = u — (A* — R). The equations resulting from the 
Lie-Poisson bracket ( p5D of the Appendix in these relative variables keep their forms and the 
condition = p is still preserved. We conclude with the following three remarks. 

Superfluid Coriolis force couples to the vortex line velocity. The Hamiltonian 
evolution equation for the superfluid velocity in the rotating frame is expressed as 

Qt< + (v: ■ V)v: = - V(p - ^|v: - v^p - ^IRH + (v; - v:) x a;* + v; X m . (41) 
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The last term is the Coriohs force and it involves the relative vortex line velocity. The curl 
of this equation yields 

dtiuj* + 2n) = curl (v; X (iu* + 2fi)) . (42) 

The form of the vortex dynamics equation is invariant under passing to a steadily rotating 
frame, and the superfiuid Coriolis force contains the renormalized vortex line velocity, rather 
than the superfiuid velocity. Therefore, this is not merely a kinematic force. The vortex 
line velocity appearing in the superfiuid Coriolis force includes the interaction between the 
vortex lines and the superfiuid component. It also includes the interaction with the normal 
component, since depends on the relative momentum density and contains finite temper- 
ature effects. The superfiuid Coriolis force is essential in the spin up problem in He-II, see, 
e.g., Reissenegger [1993]. 



Superfiuid Taylor-Proudman theorem. For steady, or slow motions and rapid rotation 
we have 

= curl (v; X 2n) . (43) 
If the rotation is uniform (Vfl = 0) and oriented vertically {fl = this becomes 

= 2\n\ (9.v;-zdivv,*) = 2|r]| id.v*,,,d.v*,y,-d.,vl,-dyv*,f, (44) 

where ( )^ denotes transpose of a row vector into a column vector. Thus, for steady, or slow 
motions and rapid uniform rotation, we find that vortex line motion is columnar. That 
is, the lateral vortex line velocity is nondivergent and independent of the axial coordinate, 
and the axial velocity decouples from the lateral motion. Therefore, under sufficiently rapid 
rotation, the superfiuid vortex filaments will straighten and become parallel to the axis of 
rotation as they approach a steady state. However, they may still undergo nondivergent 
motion in the lateral plane. This superfiuid Taylor-Proudman theorem explains why steady 
superfiuid vortices tend to be aligned with the rotation axis under rapid uniform rotation. 
The same conclusion applies, if the velocity v^ in the Hamiltonian formulation is replaced 
by the phenomenological relative velocity v^ = v^ — R. Similar considerations are discussed 
in Sonin [1987] from a more microscopic viewpoint. 

Relative total momentum is not conserved for rotating compressible flows. Since 
the Hamiltonian depends explicitly on spatial position, instead of conserving relative total 
momentum, we have the balance 



dh 



d,j: + d.rr = - ;^ =^9,|R|^ (45) 



Qx^ explicit 

where h is the Hamiltonian density in equation (^Ol). This relative momentum balance is 
the effect of centrifugal force. Here we have dropped terms proportional to p — n, since 
p = n is still preserved in a rotating frame. Consequently, the stress tensor in the relative 
momentum equation also keeps its form in passing to a rotating frame, although the total 
relative momentum is no longer conserved if the flow is compressible. 
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Appendix: Lie-Poisson Hamiltonian formulation 

Conservation of the number of quantum vortices moving through superfluid ^He (and across 
the streamlines of the normal fluid component) is expressed by 



d 
di 



JJ^u-hdS = 0, (46) 



where the superfluid vorticity cu is the areal density of vortices and n is the unit vector 
normal to the surface S whose boundary dS moves with the vortex line velocity v^. When 
CO — curl Vs this is equivalent to a vortex Kelvin theorem 

Vs-dx = 0, (47) 

which in turn implies the fundamental relation 

dtVs -veXcu^V/i. (48) 

The superfluid velocity naturally sphts into = u — A, where u = V0 and (minus) the 
curl of A yields the superfluid vorticity cu. The phase (f) is then a regular function without 
singularities. This splitting will reveal that the Hamiltonian dynamics of superfluid ''He with 
vortices may be expressed as an invariant subsystem of a larger Hamiltonian system in which 
u and A have independent evolution equations. 

We begin by defining a phase frequency in the normal velocity frame as 

dt(l) + Vn-V(f)^iy. (49) 

The mass density p and the phase (j) are canonically conjugate in the Hamiltonian formula- 
tion. Therefore, one may set u — — 5h/5p for a given Hamiltonian h and the phase gradient 
u — Vcf) satisfies 

dtu + v„ • Vu + (Vv„)^ • u = - V— , (50) 

dp 
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where ( denotes transpose, so that (Vv„)-^ ■ u = UjVvl. The mass density p satisfies the 
dual equation 



dtp + V ■ (pv„ 



„ 6h 



(51) 



Perhaps not surprisingly, the rotational and potential components of the superfiuid velocity 
Vs = u — A satisfy similar equations, but the rotational component is advected by the vortex 
line velocity vg, instead of the normal velocity v„. Absorbing all gradients into u yields 



dtA + X 00 = , 



Taking the difference of the equations for u and A then recovers equation 



as 



V(v„.u + -) 



with 



u 



A 



(52) 



(53) 



in which one uses regularity of the phase to set curlu = 0. It remains to determine from 
the Hamiltonian formulation. Including the additional degree of freedom A associated with 
the vortex lines allows them to move relative to both the normal and super components of the 
fluid, and thereby introduces an additional reactive force without introducing any additional 
inertia. This Hamiltonian approach thus yields renormalized HVBK equations. 



Proposition: Upon splitting the superfiuid velocity into = u — A (with u = V0 so that 
uj = — curl A j the (renormalized) HVBK equations in the Galilean frame of the normal fluid 
form an invariant subsystem of a Lie-Poisson Hamiltonian system, 



dl 

dt 



{fjh} with 
and Lie-Poisson bracket given by 

{f,h} = - 



f,he{M,p,S,u,A,n) 



(54) 



\6Mi 



„ r, , ^ . 5h n 5h ry^ ,^ .Sh 

{Mkd, + dkMj)-— + pdj— + Sdj— + {dkUj -Uk,j) — 

OlVIk op Oo OUk 



+ 



^dkp + I^Sfc S + ^{ukdj + Uj^k] 
6p 6S 6uj 



6h 



+ 



Ak,j 6h 



6p 
5f^ Sh 

Ok 

6n 6Ak 



6f ^ 6h 6f ^ Sh 
-^Ok- h -^Oj — 



5uk 



6uj 6p. 



(55) 



Remarks: Here M is the total momentum density, the total mass density is p and the 
entropy density is S. We shall interpret the density n later, after we develop the Hamiltonian 
equations of motion. It shall emerge that n = p is an invariant condition and, hence. 



M-nA 



(56) 



for n = p, where p = J — pv„ = Ps(vs — v„) is the relative momentum density of the 
superfluid in the Galilean frame of the normal fluid. The momentum density associated 
with the vortex fluid will be N = —nA. The Hamiltonian will be the energy E in (|^). 



Renormalized HVBK dynamics for He-II 



D. D. Holm 



14 



The Lie-Poisson bracket in the Proposition appeared first in Holm and Kupershmidt 

[1987] in a study of various approximate equations for the dynamics of multicomponcnt su- 
perfiuids with charged condensates. The mathematical structure of this Lie-Poisson bracket 
and its association with the dual of a certain Lie algebra is discussed in Holm and Kuper- 
shmidt [1987]. Our re- interpretation of this Poisson bracket introduced and studied earher 
shall now yield a extension of the HVBK equations that enables the vortex line velocity 
V( and hence the vortex reactive force and mutual friction force to be expressed at finite 
temperature. Identifying this Poisson bracket as being dual to a Lie algebra establishes that 
it satisfies the Jacobi identity, eijk{fi , {fj , fk}} = 0. The term in the Poisson bracket re- 
sponsible for the reactive vortex force will turn out to be {Ai , Aj} ^ 0. The Poisson bracket 
{vsi , Vsj} would vanish (as does {ui , Uj} — 0) and thus the reactive vortex force would be 
absent, in any Hamiltonian formulation for which A and n were not independent degrees 
of freedom from M, p, S. Volovik and Dotsenko [1979, 1980] obtain a different result and 
provide no Lie-algebraic justification. 

A Lagrangian formulation of these equations is also available. However, it involves an 
equation for 51 /Sv about which nothing is known physically. 



Corollary t^I: The Lie-Poisson bracket is equivalent to the following separate Hamiltonian 
matrix forms for the dynamical equations 









' Mjdi + djMi 




pdi 


djUi — Uj 


d 


S 




d,S 











m 


P 




djP 








dj 








_ Ujdi + Uij 





di 






'J ,* 



5h/SMj 
6h/6S 
dh/6p 

Sh/Suj 



and, upon defining N = 


-uA, 




d 


' ' 




di 


n 





ndi 




Sh/SNj 
Sh/dn 



These are individually expressed as 



(57) 



(58) 



dtS = {S, h} = - div{S Sh/5M) , 

dtu — {n,h} — —div{n6h/6N) , 

dtp = {p,h}^-diy{p5h/5M + 5h/5u), 

dtU = {u, /i} = - V( 5h/5p + {5h/5M) ■ u) + {5h/5M) x curl u , 

dt{N/n) = {(N/n), /i} = - V( 5h/5n + {5h/5N) ■ (N/n)) + {5h/5N) x curl (N/n) , 

dt(Mj + Nj) = {Mj + Nj,h}^-dkTjK 



Corollary #2: Consider a translation invariant Hamiltonian density with dependence 

h{M,p, S,n,Vs,u,A), 
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where = u — A, A = — N/n and uo = curl v^. The stress tensor T-' is expressed in terms 
of this Hamiltonian as 



where 



(59) 



^ ^ , dh dh ^dh dh ^ , , 

as m the Euler relation for pressure. 

Remark. One notes many parallels and correspondences among these equations. Note 
especially the expected similarities in the equations for u and N/n. Recall that A = — N/n, 
so that the superfluid velocity is given byv^ = u — A = u + N/n. The evolution of the 
superfluid velocity is consistently composed as the sum of these two separate dynamical 
pieces. 

Proof of the Proposition: The following Hamiltonian h (and conserved energy) will yield 
the HVBK equations in the frame of the normal fluid upon using this Lie-Poisson bracket 

h = ( d^x\ - \pvl + (M - pA) ■ v„ + eo(p, •S', - v„, u) 



(61) 

The variational derivatives of the Hamiltonian h are computed in this reference frame by 
using the thermodynamic first law (0). Namely, 

6h = J d^x[{n-^vl-A-Vn)Sp + T6S + \n-6M + {p + cm\\)-6u 

- (p + curl A + pv„ ) ■ 5A + (M - p - pv„ - pA) ■ 5v„ . 

Here we have used the velocity split ^v^ = (5u — 6A and assumed the boundary condition 
h ■ u X A = when integrating by parts. This boundary condition is satisfied identically, 
since A = A ci) in the HVBK theory. Upon substituting these variational derivatives into the 
Lie-Poisson bracket. Corollary 7^1 yields the following equations expressed in the normal 
fluid reference frame, 

dtS = {5,M = -div(5v„), 

dtn = {n, /;,} = — div(pv„ + p + curl A) , 

dtp = {p, /i} = -div(pv„ + p + curlA), 

( Hence, the condition n = p is preserved.) 

dtu = {u, /i} = -V(/i- + v„ ■ Vs) + v„ X curlu, 

( Hence, curl u = is preserved.) 
dtA = {A, /i} = n~"^(pv„ + p + curl A) X curl A , 

( Hence, V£ = v„ + p^^(p + curl A) when n = p is used.) 
dt{Mj - nAj) = {Mj - nAj, h} = - dk T/ 
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RemEirks: 

(1.) Preservation of the condition n — p by these equations allows the introduction of the 
momentum-carrying field A as an independent degree of freedom without introducing 
additional material inertia, provided the dynamically preserved condition n = p holds 
initially. This is reminiscent of the preservation of Gauss's Law by the continuity 
equation for mass conservation in a fluid plasma. 

(2.) The curl of the dynamical equation for the field A implies the vortex line velocity 

1 Sh 

Vi — — - = V + p~^curl A , where v = v„ + p~^p = J . (62) 

n oA 

The velocity v is the mass averaged velocity. The vortex slip velocity s corresponding 
to the vortex line velocity is the basis for the phenomenological reactive and mutual 
friction forces f and Rayleigh dissipation function R in the HVBK system. Namely, 

s = Vi — Vn = p~^{p + curlA) , with A = A . (63) 

As expected, this expression agrees with BK [1961] at zero temperature. Note that 
the renormalized HVBK equations introduce no new parameters. 

(3.) The corresponding equation for = u — A is then obtained as 

dt^s = - V(/i - ^ + v„ ■ v^) + V£ X , where = v + p~^curl A u . (64) 

This may be expressed equivalently in manifestly Galilean invariant form as 

dtv, + (v, • V)v, = - V( // - ^|v, - v„H + f ' , where f ' = (v^ - v,) x a; . (65) 

The term f is the Hamiltonian contribution to the reactive vortex force. This 
contribution would vanish if the vortex lines moved with the superfluid velocity. 

(4.) The stress tensor — total momentum conservation is given by summing 

Tt/ = ( Ps Vsj V^^ + pn Vnj V^) + P , t/ = Oi €klm{VsjXm) - Xj UJ^ + UJ ■ . (66) 

The divergence of defines the vortex line tension T as 

a^r/ = -a; • VA + V(cu ■ A) = -p,T + V(cu ■ A) , (67) 
In the stress tensor the pressure P is defined by the Euler relation, 

P^-eo + pp + TS, (68) 
so that in the normal-fiuid frame the pressure satisfies 

dP = pdp + SdT — p ■ d{vs — v„) — \ ■ duj . (69) 

The stress tensor Tj' = vrj^ + rj' may be derived by using Corollary ^2 for the Hamil- 
tonian formulation. 
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Implications of the HVBK vortex dynamics. The new Hamihonian formulation of 
the renormalized HVBK equations presented in the Proposition provides a formula for the 
slip velocity of a vortex line element in a turbulent superfiuid at Gnite temperature. Namely, 
for the Hamiltonian h in equation (|6TD , one finds 

s = V£ — v„ = (p + curl A) . (70) 

This formula for V£ recovers the HVBK expression in BK [1961] at zero temperature. Oth- 
erwise, it provides an extension to finite temperature of the HVBK vortex force 

f = (vl - Vs) X w , with vl = v„ - ps{as + Plj X s) , (71) 

where the renormalized vortex slip velocity is given by 

s = Vf — v„ = — (v° — v„) , for constant ps ■ (72) 
P 

The corresponding heating rate R is given by 

R = {3 — p^n + curl X) ■ u X (v^ — v„) = pps(3 uj\s x , (73) 
which is positive, as it must be. 
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